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The Modern Description of Semileptonic Meson Form Factors 

Richard J. Hill 

Fermi National Accelerator Laboratory 
P.O. Box 500, Batavia, IL 60510, U.S.A. 



I describe recent advances in our understanding of the hadronic form factors governing semileptonic meson 
transitions. The resulting framework provides a systematic approach to the experimental data, as a means of 
extracting precision observables, testing nonperturbative field theory methods, and probing a poorly understood 
limit of QCD. 
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1. Introduction: into the meson 



Semileptonic transitions of one meson into an- 
other yield important measurements of both weak and 
strong dynamics. By comparing the experimentally 
determined decay rate to a theoretical normalization 
of the relevant hadron transition amplitude at one 
or more kinematic points, elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix are determined. 
Independent of the overall normalization, the shape 
of the semileptonic spectrum provides a quantitative 
probe of underlying hadron dynamics. 

Of the six CKM elements that can be probed di- 
rectly using stable hadrons, determinations from ex- 
clusive semileptonic transitions are the most precise 
(iKisI); competitive with other determinations (IVubl, 
\Vcb\ from inclusive semileptonic B decays; \Vcd\ from 
deep-inelastic neutrino scattering; \Vcs\ from charm- 
tagged W decays); or complementary to existing de- 
terminations (I Kid I, from nuclear beta decay) lij. In 
all cases, the theoretical normalization gives a domi- 
nant error. The experimentally determined spectrum 
can be used both to test the nonperturbative methods 
used in determining this normalization, and to opti- 
mize the merger of theory with experiment. 

For the study of hadron dynamics, the use of a vir- 
tual W boson probe in semileptonic transitions can be 
viewed in analogy with the use of a virtual photon in 
deep-inelastic scattering. The "known" and theoret- 
ically "clean" weak (electromagnetic) physics is used 
to probe the "unknown" and theoretically "messy" 
strong physics of the meson (proton). For the case 
of semileptonic transitions, in addition to the virtu- 
ality of the exchanged boson, we can envision "dial- 
ing the knobs" of the initial- and final-state meson 
masses. To be precise, consider the situation pictured 
in Fig. ^ A Qq pair with quantum numbers of the 
"heavy" meson H is created at spacetime point x, in- 
teracts with a flavor-changing current at y, and the 
q'q pair with quantum numbers of the "light" meson 
L is annihilated at z. By suitable manipulations ^, 



we can extract hadronic transition amplitudes (form 
factors) from this correlation function. The entire pro- 
cess takes place in the complicated QCD background 
with soft and hard gluons interacting with everything, 
qq pairs popping out of the vacuum, etc. On this 
fixed background, we can compute a correlation func- 
tion with any quark masses we desire; in fact Nature 
has chosen a few fixed values (and presumably ob- 
tained the correct answer), and we shall be content 
with these. ^ 



Q 




P 



Figure 1: Correlation function from which the 
semileptonic form factors are extracted. Q, q' and q are 
quark labels, p'', and q^ are momenta. 



In this way, a rather complete exploration of the 
"1-body" semileptonic topology, summarized by in- 
variant form factors F(mH,mL,q^), is possible. A 
similar analysis can be applied to "0-body" leptonic 
decays, and to "2-body" hadronic decays. The for- 
mer case is rather simple: by kinematics there is only 
the single "knob" of the meson mass, and the result 
is summarized by a single number, the decay con- 
stant /(m^f). The latter case is rather complicated, 
with many different topologies contributing to a typ- 
ical physical process • 

The remainder of the talk is organized as follows. 
Section [3 reviews the first-principles knowledge we 



Fourier transform, x 



p, z 



p', and tune p 



^(Unquenched) lattice simulations study essentially the same 
object. However, in contrast to the fixed background described 
here, the "valence" quark masses injected by the currents are 
coupled to the dynamical "sea" quark masses. When extrapo- 
lated to the physical masses, the results must of course be the 
same. 
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Table I Maximum |2;(t, to)| throughout semileptonic 
range with symmetrizing choice to = t+i^ — -^/l — t-/t+) . 
Process CKM element U max 







Vud 


3.5 X 10 


B - 


^ D 


V,b 


0.032 


K - 




Vus 


0.047 


D - 


^ K 


V,s 


0.051 


D - 


-* n 


Vcd 


0.17 


B - 


-* TV 


V^t 


0.28 



Figure 2: Mapping ^ of the cut t plane onto the unit 
circle. The semileptonic region is represented by the blue 
line. 



have about the form factors, following just from kine- 
matics without dynamics. Pseudoscalar-pseudoscalar 
transitions between "heavy-light" , nonsinglet mesons 
are particularly simple and are the main focus. 
Rigorous power-counting arguments provide the basis 
for a powerful expansion based on analyticity. Sec- 
tional illustrates how the experimental data is simpli- 
fied by making use of this expansion. In particular, 
we find the remarkable conclusion that in terms of 
standard variables, no semileptonic meson form fac- 
tor has ever been observed to deviate from a straight 
line. Given that the form factors are indistinguishable 
from straight lines, if the shape of the semileptonic 
spectrum is to provide insight on QCD, it must be 
through the slope of the form factor; in fact, a clear 
but unsolved question in QCD translates directly into 
the numerical value of this slope in an appropriate 
limit, as described in Section 0] Phenomenological 
implications in the _B — > tt system are considered in 
Section |S1 The methodology described here provides 
a convenient framework in which to understand pre- 
cisely what measurements in the charm system can, 
and cannot, say that is relevant to the bottom sys- 
tem, as discussed in Section Section [7| outlines the 
extension to pseudoscalar- vector transitions. 



2. Analyticity and crossing symmetry 



An oft-cited downside of old and well-known 
dispersion-relation arguments is that the results are 
too general, and do not make specific predictions for 
detailed dynamics. In fact, precisely these properties 
make them useful to the problem at hand — it is essen- 
tial to make some statement on the possible functional 
form of the form factors, yet we do not want to make 
assumptions, explicit or implicit, on the dynamics. 

The analytic structure of the form factors can be 



investigated by standard means. Let us focus on 
the form factors for pseudoscalar-pseudoscalar transi- 
tions, defined by the matrix element of the relevant 
weak vector current, {q = p — p') 

{Lip')\V^\H{p)) 



F+iq'){p^'+P 



'■H 



'Foiq^ 



TO? „ 

— -q'' . 



To ensure that there is no singularity at q^ 
form factors obey the constraint 

F+(0) = ^^o(0). 



(1) 

0, the 

(2) 



Ignoring possible complications from anomalous 
thresholds or subthreshold resonances, to be discussed 
below, the form factors F{t = q^) can be extended 
to analytic functions throughout the complex t plane, 
except for a branch cut along the positive real axis, 
starting at the point t = t+ \t± = (m^f ±TOi)^] corre- 
sponding to the threshold for production of real HL 
pairs in the crossed channel. By a standard transfor- 
mation, as illustrated in Figure |21 the cut t plane is 
mapped onto the unit circle \z\ < 1, 



z(t, to) 



(3) 



where to is the point mapping onto z — 0. The iso- 
lation of the semileptonic region from singularities in 
the t plane implies that \z\ < 1 throughout this re- 
gion. Choosing to = ^+(1 ^ -y/l — t-/t+) minimizes 
the maximum value of \z\; for typical decays these 
maximum values are given in Table D 

Since the form factor is analytic, it may be ex- 
panded, 



m = 



1 



P(t)0(t,to) 



'^ak{to)z(t,to)'' , 



(4) 



fc=0 



^The nonsinglet restriction ensures that only a single topol- 
ogy is relevant as in Figure [Tl 



*For a general discussion, see e.g. 01 . For earl; 
applications to semileptonic form factors, see 



work_on 
IC 
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where P and cj) will be explained shortly. From Ta- 
ble it is apparent that if some control over the coef- 
ficients flfe can be established, the expansion is rapidly 
convergent. 

To bound the coefficients appearing in Q), we con- 
sider the norm, 



OO 

fc=0 

°° dt 



Til J Z 



1 



TT 



+ 



t - t. 



-\P^F\\ (5) 



By crossing symmetry, the norm can be evaluated us- 
ing the form factors for the related process of HL 
production. 



2.1. Subthreshold poles, (absence of) 
anomalous thresholds, and a choice of P 

For some hadronic processes, it may happen that 
subthreshold resonances occur in the production am- 
plitude, which must be properly taken into account. 
Particles lying below threshold are hadronically sta- 
ble, so that ignoring higher-order weak and electro- 
magnetic corrections, they are described by simple 
poles. The canonical example is the B* pole appear- 
ing in the vector channel for _B — s- tt. Such poles could 
in principle be simply subtracted, but doing so re- 
quires knowledge of the relevant coupling appearing 
as the coefficient of l/(m^, — t) in the dispersive rep- 
resentation of the form factor. Armed with only the 
knowledge of the pole position, this pole can instead 
be removed by multiplying with a function P(t) with 
a simple zero at m^, . Requiring also that the func- 
tion satisfy \P\ 
phase. 



1 along the cut, up to an arbitrary 



-IF 



(6) 



with z as in (|^. 

It may happen in some processes that "anoma- 
lous thresholds" appear. The relevant aspects of this 
technical subject can be summarized as follows: an 
anomalous threshold can occur in the spacelike region, 
t < 0, only if H and L are unstable, i.e., rnn or l > 
mx + my for some hadrons X and Y coupling to H 
or L. Similarly, an anomalous threshold can occur in 
the timelike region, t_ < t < t^, only if rn-j^^^^ > 
m\ + niy- ^ This explains the priveleged position of 



ground-state heavy-light pseudoscalar mesons. The 
ground state meson with given flavor quantum num- 
bers is necessarily pseudoscalar .11.] . The mesons H 
and L are therefore the lightest hadrons containing 
their respective "heavy" quarks, so that in particular 
iTT-H or L < 'rnx + rny and or l < "^x + ^^r any 
hadrons X and Y containing the same heavy quark. ^ 
(For the present purposes, this can be taken as the 
definition of a "heavy-light" meson.) Fortunately, and 
for related reasons, these mesons are easily produced 
and studied experimentally, and in lattice simulations. 



2.2. Unitarity and a choice of 

Nothing in 101) or (jsj yet singles out a choice of 0; 
indeed any analytic function will work, (of which there 
are many!). A default choice is determined from argu- 
ments based on unitarity: by an appropriate choice of 
(p, the norm can be identified as a partial rate for some 
inclusive process that is perturbatively calculable. In 
particular, from 

n'^'iq) = i J d^xe"^ '={0\T{V^{x)V''^{0)}\0) 

= {q^q" ~ 5'"'g')ni(g2) + q^q'^^^i^q^) , (7) 

unsubtracted dispersion relations can be written for 
the quantities (Q^ — —q^) 



Noticing that for t > t+, (ry an isospin factor) 

[(i-<+)(t-t_)p/2 



(8) 



487r 



r^t+t^ [(f-t+)(f-t_)]V2 



167r 



F+{t)\^ < Imni(t), 
\Fo{t)\^ < Imno(i),(9) 



shows that an upper bound on the norm can be es- 
tablished by choosing [recall that |z| = 1 along the 
integration contour in ©] 



V t+-t ( z{t,0) 



487r {t+ - to)^/* V -t 



^In geometrical language, this can be related to the state- 
ment that a triangle cannot have more than one obtuse an- 
gle 01 . Strictly speaking, for "heavy-to-heavy" transitions such 
as D — > X and B —> D, additional Zweig-suppressed topologies 
can also lead to anomalous thresholds, related to processes such 



as Ds~^ — > Di^Tr". As indicated by the small branching frac- 
tion, such effects are highly suppressed; for further discussion 
and references, see Section 171 

^The same is not true for "heavy-heavy" systems; e.g. a QQ 
pair has mass ~ (2mQ)^ ~ ^'"q' compared to -|- rUg ~ 
2m?j for a pair of (Qq) mesons. 



fpcpOe 324 



4 



Flavor Physics and CP Violation Conference, Vancouver, 2006 



t- - 1 



4'Fo{t,to) = 



r]t+t^ y/t+ - t ( z{t,0) 



16n (t+ -to)i/4 \^ -t 



z{t,-Q^)\ (z{tM) 



-Q2 -t J \tQ-t 



-1/2 



i_ - t 



(10) 



The choice of subtracted dispersion relation in JSJ 
leads to a "default" choice for (j) in fTn| . With 
this choice, power counting shows that ak/ao and 
^ki^k/^o) ^'^ ^'^^ scale as powers of large ratios such 
as Q/Apcp, or toq/Aqcd when a heavy-quark mass 
is present [53 ■ This ensures that there is no paramet- 
ric enhancement of the coefficients Ofc that could offset 
the smallness of z'^ in the series In fact, at suffi- 
ciently large k the coefficients must decrease in order 
that the sum of squares converge. These properties 
[analyticity, and Uk/ao ~ Sfc'^fc/^o ^(1) ] 
that is required from the choice of (f>. The "physical" 
prescription following from |(SJ), © and (|1U|) automat- 
ically ensures that this is the case. 

The original motivation for considering the operator 
product expansion (OPE) in (|7|) is to place a restric- 
tion on the coefficients in (O according to (x — xf+ o 
as appropriate) 



10 
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-e- 



10 

t (GeV^) 



20 



-0.2 



L_ 



-z 



0.2 



k=0 



(11) 



However, in order that an OPE expansion for x(Q^) 
converge, (or TOq -I- when a heavy-quark is 
present) must necessarily be large compared to Aqcd • 
This results in a bound that is typically overestimated 
by some power of the large ratio of perturbative to 
hadronic scales. In practice, the numerical value for 
the bound itself is largely irrelevant. What is im- 
portant is that the choice of (j) which it motivates has 
the desired properties. 

Having chosen a "default" 0, for definiteness, we 
will also take Q = in (|10l) as the "default" choice, 
and where a particular choice is necessary, to = t+i^~ 

Vi-i-A+)- 



3. What the data say 

Table m can be used to predict the level of precision 
at which slope, curvature, and higher-order correc- 
tions can be resolved by the data. With the "default" 
values of Q and to, Tablellllshows the results for ai / ao 
obtained from data. Except where indicated, modes 



Figure 3: Experimental data for F+ in _B — > nlv, as a 
function of t, and for the function P(f>F+ as a function of 
z. Partial branching fractions have been converted to 
values of the form factor at the midpoint of each bin in t 
from [l^. Units on the vertical axis are arbitrary 
(approximately normalized to unity at t = 0). 



related by isospin are combined. For the K n 
case, the results of 0, 0, 0, 0| were presented as a 
simple quadratic Taylor expansion of the form factor 
about t = 0. ^ These results have been converted to 
the quadratic z parameterization in Q), by identifying 
the Taylor series at t = 0, and propagating errors lin- 
early. For B ^ D, the results of [iJl were presented 
in terms of a parameterization obtained by expand- 
ing and P as a Taylor series in z |^. The result in 
Table llll is obtained by converting to the linear z pa- 
rameterization in with three subthreshold "i?*" 
poles located at m = 6.337, 6.899, 7.012 GeV Hil, 



^It is desirable to fit the data dire ctly to jl), to avoid biases 
introduced by the truncated t series l23l . 
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Table 11 Linear expansion coefficient ai/ao from 10} at 
to = t+{l - ^1 - t^/t+) and = 0. 



Process 


ai/ao 


Reference 


B ^ D 


-2.6 ± 2.3 


[13] 


K+ ^ tt" 


-0.2 ±0.2 


[14j 


Kl 7r± 


-0.5 ±0.2 


[15] 




0.0 ±0.3 


[16j 




-0.2 ±0.2 


[171 


D ^ K 


—2 7 ± 5 ± 4 


[181 




-2.2 ±0.4 ±0.4 


[19j 




-3.2 ±0.5 ±0.2 


[2i2j 


D ^ n 


-2.3 ± 0.7 ± 1.3 


[18j 




-1.6 ± 0.5 ± 1.0 


[20] 


B n 


-1.3 ±0.6 ±2.3 


[2JJ 




-1.9 ± 0.3 ± 1.1 


[12j 




-1.3 ±0.8 ±2.2 


[221 



and then identif ying the co efRcients in a Taylor series 
at t = t+{l- ^1 - The results for D ^ K, 

— > TT and B —t t: were obtained by fitting the linear 
z parameterization in 10} to the data. A second error 
is included by redoing the fits with the quadratic z 
parameterization, subject to the conservative bound 
|a/c/ao| < 10. Due to the smallness of z for pion beta 
decay, 7r+ tt" (cf. Table P), the slope in this case is 
orders of magnitude from being measured experimen- 
tally m. 

At least through linear order, there is no evidence 
of anomalously large coefficients that could upset the 
power counting. While it would be desirable to push 
to the next order and examine the size of 02/00, com- 
parison to data establishes the remarkable conclusion 
that form-factor curvature has not yet been seen in 
any semileptonic transition. In fact, for many cases, 
a form factor slope has yet to be measured. An ex- 
ample of the transformed form factor is illustrated for 
B TT in Fig. 

From the amusing coincidence that z^^ j z^'^ ^ 
\Vcd\l\VcA. and z^^lz^^ ^ it turns 

out that the higher statistics of the Cabibbo-allowed 
modes [B ^ D, D ^ K) arc offset at linear order 
by the smallness of z. It is thus likely that curva- 
ture will eventually be measured first in the Cabibbo- 
suppressed modes (i? — s- tt, D n). 

The results in Table are by no means the fi- 
nal word on these quantities, but illustrate the main 
point, that there is no sign that the z expansion is 
breaking down. It is also easy to see that unitarity has 
very little impact. For example, for i? ^ tt, the bound 
on taken from the OPE at Q = is overestimated 
by a factor ^ (?Tif,/AQCD)^ _53J. Taking for definite- 
ness, F+{to = 16GeV) « 0.8, the unitarity bound 
tells us only that Y.k "fc/oo - 2500 d IM El- For 



i? ^ D, at Q = with the approximate symmetry re- 
lation i^_^(t_) w [niB+'mjj) /2^mBTn£), and including 
three subthreshold poles as in © , the unitarity bound 
is overestimated by a factor ^ {mi,/mc)^ and yields 
Efc flfc/ao ~ 9000 11 0- While these bounds can be 
improved somewhat by subtracting off subthreshold 
poles, extending isospin SU{2) to SU{?>) flavor sym- 
metry, or by lowering Q^, all of these modifications 
introduce their own uncertainties. ® 



4. A fundamental question 

Given that the form factors (after extracting stan- 
dard kinematic factors, and expressing them in terms 
of the appropriate standard variable) are so far indis- 
tinguishable from a straight line, it is apparent that 
any insight to be gained from the shape of the form 
factors, whether it be tests of nonperturbative meth- 
ods, inputs to other processes, or more fundamental 
questions about QCD, must be based in first approx- 
imation on the slope of the form factor. In fact, this 
quantity does provide a clear test of lattice QCD, is an 
important input to hadronic B decays, and in an ap- 
propriate limit can provide the answer to a longstand- 
ing open question about the QCD dynamics governing 
form factors. 

It is convenient to define the physical shape observ- 
ables in terms of the form factor slopes at t = , 



1 



5 = 1- 



-(0) 



i dF^ 
dt 

9 

mi 



F+(0) 
F+(0) + F_(0) 
F+(0) 



dF+ 



dt 



_ dFo 




t=o dt 





(12) 



The quantities (3 and d depend only on the masses 
of the mesons involved. ^ Being physical quantities, 
they are independent of any renormalization scale or 
scheme. As discussed in the introduction, these quan- 
tities take definite values for all mn and mi,, values 
that are accessible experimentally at the fixed masses 
ttItt, ru]^, mo and ttlb- 



"For modes such as -B — > D, the incredible smallness of z, 
and the judicious use of heavy-quark symmetry, allows even 
very conservative unitarity bounds to guarantee few-percent 
level accuracy by keeping only the linear term in BB- 

^Recall that we consider mesons with a fixed light spectator 
quark, which for simplicity in the discussion is assumed mass- 
less. The meson mass is therefore in one-to-one correspondence 
with the heavy (non-spectator) quark mass. 

l°For some studies of 5 - 1 = F_(0)/F+(0), also called §(0), 
in the early literature of light-meson form factors, see the re- 
view l2Hll . T he posi tive sign for ^(0) predicted in a number of 
models, e.g. I3(ll3ll . is in disagreement with current data. For 
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Due to the kinematic constraint (0), the difference 
of form factor slopes is particularly simple. Firstly, 
for rnu = m^, 



(13) 



This is the statement of current conservation (F_ — 0) 
in There are three further distinct limits that we 
can consider: mH,mL > Aqcd, mH,mL < Aqcd, 
and niH S> Aqcd 3> mi^. Each limit provides valu- 
able insight, and we consider each in turn. 



4.1. mH,mL > Aqcd: HQET 



CO 



-1 



2 4 

( GeV ) 



Figure 4: 5 as a function of light-meson mass for fixed 
heavy-meson mass mn = ms. The shaded band gives 
the prediction from HQET valid for mn 1 GeV. The 
dot with error bars is from experimental data as 
discussed in the text. 

For mH,rnL 3> Aqcd, it is convenient to express 
the form factors F± in terms of reduced amplitudes, 
h± where the dominant heavy-quark mass dependence 
is extracted: [cf. (^] 

F+{p + p')+F^{p-p') 

= y/mnmL [h+{v + v') + h^{v - v')] , (14) 

where the meson velocities are mnv = p, m^v' = p' . 
With these definitions, before any approximation. 



5{mL,mH) 



2niL 



h__ 

h+ 



(15) 



rriff+mi h+ 



some early work on heavy-to-light meson form fac tors, see |32| . 
The prediction limmg^txi S{m,Tr ^mg) = 2 of |32| also appears 
difficult to reconcile with present data, see below. 



At leadin; 
function 



5{mL,mH) 



ower, h+ 
and ? 

2mL 



« 1^, the universal Isgur-Wise 
0, so that 



■[l-KC'(as,A/TOL,A/mff)] , 

rrLH + mL 

[TO//,mL>A]. (16) 

Radiative and power corrections to h± can be ana- 
lyzed systematically with heavy-quark effective field 
theory (HQET) [SJ- From (O, we see that in the 
regime where mn, 3> Aqcd, S can take any value 
between zero and unity. For fixed mn = rns, Fig. 01 
shows the allowed range of 6 in the regime of mi, 
where this expansion is applicable. Results for sub- 
leading corrections are taken from j^, with errors 
estimated by varying the renormalization scale by a 
factor of two, and assigning 100% uncertainty to the 
power corrections. The data point for _B ^ tt is from 

4.2. mH^rriL < Aqcd: CHPT 



to 




Figure 5: 5 as a function of heavy-meson mass for fixed 
light meson mass rriL = rn-jr- The shaded band gives the 
prediction from CHPT valid for ruH 1 GeV. Dots with 
error bars are experimental data as discussed in the text. 

For mH,mL Aqcd, the form factors can be ex- 
panded in powers of the small ratios of masses and 
momenta relative to the QCD scale (since the energies 
involved in semileptonic transitions are bounded by 
meson masses) . In this regime of Chiral Perturbation 



Since we are considering quantities at maximum recoil, the 
large recoil parameter, {v ■ i;')max = (m^ + mj^)/2mLmH , can 
upset the power counting when mj{ — > oo at fixed mi^ . For arbi- 
trarily small mi^/rriH, the limiting value 5 = 2mj^ j (mu -^rai,^ 
is obtained in the limit where mi^ — > cxd, ran oo, mi^/rriH = 
constant . 
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Theory (CHPT) [33, |36| , 1 - (5 is of order / Kf^^^. 
The leading contribution is given by 

2{m% - ml) 



= 1 - 



[Lli„)~2Llif,) + ...], 
[TO^/,mi<A], (17) 



where the ellipsis denotes a known kinematic function 
which cancels the renormalization-scale dependence of 
the low-energy constants g{fj.). Here Fq (x is re- 
lated to the pion decay constant. The deviation of 5 
from unity is predicted by the sign of the combina- 
tion Lg — 2L^, which is empirically found to be pos- 
itive. The band in Fig [3 shows the allowed range of 
S, where for illustration we take |[3Q] Fq = 0.088 GeV, 
10'^L5(to,;) = 2.2 ± 0.5 (determined from the ratio 
Ik/ f-w) and 10^Lg(m,,) = 7.4 ± 0.7 (determined from 
the electromagnetic charge radius of the pion). The 
value of 5 for 7r+ is taken from El- The 

data point for D t: la discussed in Section El and 
the data point for S ^ tt is the same as in Fig. 0] 



4.3. mj/ > A > niL. SCET 



can thus be analyzed using a simultaneous expansion 
in K/ruH and A/E. The soft-collinear effective the- 
ory (SCET) framework has been developed to study 
this regime 

m E El- The leading descrip- 
tion of the form factors for pseudoscalar-pseudoscalar 
transitions is, up to corrections of order as{mH) and 
K/mH k/E, (2iyiil0|4lilj43 

F+{E) = ^[C{E)^(^^-\^H{E)\, 



^Fo{E) = V^[C{E)+H{E)], 



(18) 



where it is more natural to work here in terms of the 
light-meson energy, related to the invariant momen- 
tum transfer t in ^ hy t = m\ + m\ — 2m^f£'. 
The construction of SCET is more intricate than ei- 
ther HQET or CHPT, due to the nonfactorization 
of large and small momentum modes in some pro- 
cesses Ji^^^. Apart from scaling violations related 
to this phenomenon, the functions C and H both have 
an energy dependence C, ^ H ^ E^^ . 

From expression (|18|l it is straightforward to see 
that, independent of any model assumptions, [23 | 




m, ( GeV ) 



Figure 6: Contours of constant S in the rriL — mn plane. 

The HQET description in Section im breaks down 
when the light meson becomes light (Fig. Simi- 
larly, the CHPT in Section breaks down when the 
heavy meson becomes heavy (Fig. . 

In this regime, since we are concerned with the 
point at maximum recoil, the light-meson energy nec- 
essarily satisfies E ~ m///2 Aqcd- Observables 



^^The slope parameters there are related to those in 1121 by 
1 + 1/(3-5^ \+{m]^ - ml)/ml , 1//3 Ao{m|^ - -mD/ml. 



1 _ dln{C + H) 
dlnE 



l + 0{a,,A/mH) 



(19) 

From the 1/E'^ asymptotic behavior of C and H, l/f3 
should approach unity as — > oo at fixed ^ 
Aqcd- Similarly, 



5{mL,mH) 



2H 



C + H 



0{as,A/mH) 



E=mH/2 



[niH > A > uil] 



(20) 



Thus the question of the asymptotic limit, 
lim„i^^oo i5(r7ii, TOjj), is the same as the ques- 
tion of which contribution, the "hard" function H, or 
the "soft" function C, dominates in this limit. 

The results summarized by (|16f) . H17|l and H20|) . can 
be unified as in Fig. which displays three equal-^ 
contours in the HQET and CHPT regimes. The obvi- 
ous question is: what happens to these contours as we 
pass outside of these regimes? We know that S can- 
not depend strongly on m^, since this dependence is 
due solely to scaling violations and power corrections. 
A smooth chiral limit of the form factors also implies 
that S cannot depend strongly on m^. To further con- 
strain these contours, experimental data can be used 
as a quantitative probe in this regime. 



^For explorations along different lines, see 
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5. Implications for phenomenology 

We focus attention on the prototypical heavy-light 
process, B iiiv. The methodology described in 
Section El can be used to extract as much information 
from experimental and lattice data as possible, with 
believable error estimates. 



estimates used to bound other theoretical errors en- 
tering a \Vub\ determination — e.g. power corrections, 
pcrturbative matching corrections, or discretization 
errors in lattice calculations of the form factor. Quan- 
titative investigations such as in Table|n|and SectionEl 
give us confidence that as far as the bounds are con- 
cerned, "order unity" really means order unity. 



5.1. Precision measurements: K 



ub 



5.2. Inputs to hadronic B decays 



375 




20% 
15% 
10% 



-5% 
-10% 
-15% 
-20% 



5 10 15 20 25 

t (GeV^) 

Figure 7: Experimental error on |Kji)|, with theoretical 
input of the form factor normalization at a given value of 

from [5^ . The plot assumes that the fit yields the 
central value \Vub \ = 3.7 x 10~^. 

Knowing that the true form factor is given by one 
of the restricted class of curves in Q allows for maxi- 
mal use of the available experimental ^'^ and lattice 
data. Figure \7\ shows the minimum error obtainable 
for I Vub I using present data 0, , with a 

form factor determination at a given value of t [s^ . 
The dark, medium and light bands correspond to in- 
creasing levels of conservatism for the size of the co- 
efficients appearing in Q): J2k'^k ^ O.Of,0.f,l, nor- 
malized relative to the default unitarity bound. At 
the point to = 16GeV^, these values correspond to 
Efe«fc/ao - 25,250,2500. 

As Fig. [71 illustrates, theory inputs at either very 
large or very small t are not as effective as for moder- 
ate values, say t — 10 — 20 GeV^, which are within the 
range currently studied with unquenched lattice simu- 
lations |5J, [53, ISal • For extreme values of the bounds, 
e.g., allowing coefficients in the expansion I^J to be as 
large as Uk/ao ~ 50, the error begins to increase, as 
the lighter bands in the figure show. 

For consistency, the heavy-quark power counting 
used to establish bounds on the form factor shape 
should be at least as robust and conservative as similar 



^"^For the status of B — > niu measurements see l5ll 
'^^For recent revi ews and references for lattice form factor 
determinations see l57ll5gl . 



Semileptonic decays provide a robust value for the 
form factor normalization (times |Kib|), a key input 
to factorization analyses of two-body hadronic de- 
cays HiSIIlil. From jU 



16 



10^\VubF+{0)\ = 0.92 ± 0.11 ± 0.03 . 



(21) 



A dominant uncertainty in many factorization pre- 
dictions is the normalization of the hard-scattering 
contribution to the form factor, commonly expressed 
in terms of an (inverse) moment of the B meson wave- 
function, Ar: 



5{mT,,mB) 



Nc mBXBF+{Q) 



(22) 



For example, the "default scenario" inputs of |63j give 
while the "S2 scenario" gives a central 



0.15 



+0.10 

-0.05' 



value 6 ~ 0.29. The semileptonic data can help 
pin down this number. From 53J, using experimental 
data for F+ from [Ulilllll^ determines 

B ^TT : 

l + l/P-S = 1.5 ±0.6 ±0.4, (23) 



and the lattice value /3 = 1.2 ±0.1 from |5J,|55j allows 
extraction of S. 

If 5 is monotonic as a function of m^, the analysis 
of Section 14.11 shows that 

S{mL,mB) < 0.35±0.03, [tul < 1 GeV; monotonicity] 

__(^) 

A significantly larger value of 6, e.g. S ^ 1 |65L IGo . l67| . 
would require a dramatic behavior of the extrapolated 
curve in Fig. ^ 



I'^This may be compared to the bound 63] 10^\VubF+{0)\ < 
1.22 obtained from an assumption of form-factor mono- 
tonicity and experimental data in |21| . and the value |64| 
10^|Vul,F+(0)| = 0.83 ± 0.16 obtained from fits of model pa- 
rameterizations to the same data. 

^^Values are at tree level, as(\/mi,A) 0.35, asymptotic dis- 
tribution amplitude for the pi on, and the errors shown are from 
the remaining inputs used in RiS^. 
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6. What's charm got to do with it? 



Similarly for D ^ tt decays, 



Charm decays provide a direct probe of the "in- 
teresting" regime pictured in Fig. They provide 
an important test of lattice measurements for heavy- 
to-light form factors, and a quantitative test of the 
power-counting used to bound the form factor shape 
in other processes such as B — > tt. 

6.1. Fundamental questions 



CO 



-1 



( GeV ) 



Figure 8: 5 as a function of light-meson mass. The 
upper (black) points are for fixed ttih = m_D, and the 
lower (red) points for fixed mn = m_g. 

Depending on whether the soft or hard component 
of the form factor dominates in the limit mn oo, 
the difference in form factor slopes H12|l tends to (5 
or S — > 2. The latter would require that the curves 
such as in Fig. ^ turn upward at small m^, for suffi- 
ciently large mn- If this happens, and unless a new 
large scale is dynamically generated by QCD, at which 
a turnaround in the curves would occur, some evi- 
dence for this behavior should be evident in the charm 
system. Precision measurements here directly probe 
the dangerous region illustrated in FigEl For D ^ K 
decays. 



D K : 

I + 1/(3 -5 



1.03 ±0.09 ±0.11 [ig 
0.94 ±0.07 ±0.10 [13 
1.13±0.10±0.12 113.(25) 

Here the first error is experimental using the linear 
z parameterization (just ao, oi) in The second 
error is a conservative estimate of the residual shape 
uncertainty, obtained by allowing an extra parameter 
in the fit (oq, ai and 02), with |afc/ao| ^ 10. Combin- 
ing this with the lattice value (3 — 1.8 ± 0.1 yields the 
D ^ K data point in Fig. |S| 



1 + 1//3- 



5 = 1.3 ±0.4 ±0.5 [13 
0.9 ±0.2 ±0.3 [23 



(26) 



with the errors as above. Combined with the lattice 
value (3 — 1.65 ± 0.10 yields the D ^ tt data point in 
Fig. m and Fig.il 

Our current combined knowledge of form factor 
shape from B and D decays is illustrated in Fig. |S1 
So far the data do not indicate surprises in either of 
the curves when extrapolated into the region ^ 
1 GeV. It will be interesting to probe this region with 
more precision when further data becomes available. 

6.2. Lattice, experiment and 
parameterizations 

Charm decays provide important tests of nonper- 
turbative methods used to evaluate hadronic matrix 
elements. When comparing the results of lattice QCD 
with experiment, it should be kept in mind that the 
kinematic regions that are studied with best precision 
are different for the lattice (large t) and experiment 
(small t). Also, the manner in which chiral extrap- 
olations are performed to reach physical light-quark 
masses imply that it is difficult to present lattice re- 
sults in terms of uncorrelated values of the form factor 
at different t values. In practice, the results are gen- 
erally presented in terms of a parameterized curve; 
to make a definitive comparison to experiment, it is 
essential that the chosen parameterization doesn't in- 
troduce a bias. The ideas described in Section|21allow 
a systematic approach to this problem ^|. The re- 
mainder of this subsection points out pitfalls that can 
occur with some of the simplified parameterizations 
in common use. 

The starting point for many parameterizations is a 
more pedestrian but rigorous approach to analyticity, 
which implies the dispersion relation: 



F+{0) 



1 



1 



1 — a 1 — t/mjjt 



dt 



,ln\F{t') 

t' -t 



(27) 

where a distinct mn* pole appears below threshold 
for heavy-to-light decays such as B ^ n and D ^ K 
(and almost for D — > tt). The first interesting test 
is to see whether just the ruH* pole can describe the 
data, 



vector dominance : 



1 



In fact this "vector dominance' 
itly ruled out by the data [l2.ll 



- t/m%, ■ 

(28) 

model can be explic- 

[11 [13, [13, [H [13, so 
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that inclusion of the continuum contribution in (|27|) 
is essential. 

From a dynamical point of view, in order to ob- 
tain the 1 / E"^ dependence appearing in heavy-to- light 
form factors, as in H18|l . it is necessary that the con- 
tinuum integral in 1)27(1 play a significant role. This 
can be treated in a model independent way by break- 
ing up the integral into a sum of effective poles, 
and using power-counting estimates to establish rea- 
sonable bounds on the coefficients of these effective 
poles |2^|H3- In the first approximation, the contin- 
uum integral is represented by a single effective pole, 
and two parameters are necessary to describe its loca- 
tion and strength relative to the H* pole. Since this is 
one more parameter than is easily measured from the 
data, various suggestions have been made for elimi- 
nating one of these parameters. 

The "single pole model", where both the H* pole 
and the continuum integral are represented by a single 
pole that is allowed to float, 

single pole : 



(29) 

is also ruled out by the data, although in a slightly 
less direct way. Although the curve can be made to 
fit, the pole position is forced to take an unphysical 
value, significantly below both the H* pole and the 
continuum. The "modified pole model" , 



Table III Maximum to)| throughout semileptonic 
range with symmetrizing choice to = t+i)- — -y/l — 





Process 


2 max 


D 


^ K* 


0.017 


D 


P 


0.024 


B 


D* 


0.028 


B 


P 


0.10 



6.3. Testing convergence 

The general parameterization Q provides a system- 
atic procedure for estimating how many terms should 
be resolved by data at a given level of precision. Rig- 
orous bounds are placed on the coefficients by using 
crossing symmetry to analyze the production form fac- 
tor, either via unitarity arguments, or through power 
counting of contributions from different momentum 
regions. Since the latter estimates yield constraints 
that are so much more powerful than can be safely es- 
timated by unitarity, it is important to check wherever 
possible that large "order unity" numbers don't ap- 
pear. As illustrated by Table UTI the available semilep- 
tonic data reveal no surprises. 

7. Decays of pseudoscalar to vector 
mesons 



modified pole : 

+ (l-t/m|,.)(l-«poict/™l,.) ' 

(30) 

holds a similar status. For D ^ K and D n, 
((^ and 1211) clearly show that l-f - (5 w 2 is not 
valid for charm decays, as necessary for the motivation 
for the simplification l(30() proposed in . Although 
the form 1(30(1 can be made to fit the data, there is 
no obvious physical interpretation for the resulting fit 
parameter; in particular ctpoie obtained in this way has 
no direct relation to the physical a defined in ((27(1 . 

It should be kept in mind that unless there is a 
physical meaning that can be given to the parameter 
being studied, there is no guarantee that different ex- 
perimental or lattice determinations will converge to 
any one value for this parameter. It is therefore un- 
clear what to make of discrepancies appearing when 
different datasets are forced to fit models such as (1^ 
or 1(30(1 '6^, 'T^I . The situation is especially dangerous 
for comparing lattice and experiment, since the range 
of t that is emphasized is different in the two cases. 
These pitfalls are easily avoided by working with a 
general parameterization such as 10} that is guaran- 
teed to contain the true form factor, and by comparing 
physical quantities, such as in 1(25(1 . 1(26(1 . 



Pseudoscalar-pseudoscalar transitions hold a priv- 
eleged position, from a first-principles simplicity point 
of view, from the lattice point of view, and from the 
experimental point of view. Pseudoscalar-vector tran- 
sitions, while accompanied by new complications, are 
however important backgrounds to the pseudosclar 
mode, provide alternative extractions of CKM param- 
eters, and yield important constraints on radiative 
and hadronic transitions. This section briefly outlines 
the implementation of the ideas in Section |21 to the 
pseudoscalar-vector case. 

The most obvious complication is the multiple in- 
variant form factors that accompany the vector parti- 
cle. Less obvious complications involve modifications 
to the analytic structure of the form factors due to the 
unstable nature of vector mesons, and the possible ex- 
istence of anomalous thresholds. When these occur, 
they will encroach on the gap between the semilep- 
tonic region and singularities. The effects of such 



"There have been numerous studies aimed at reproducing 
the P — > V data by means of sy mmetry argu ments or proposed 
generating resonance structures l7ll \7i l73l . The focus of the 
present talk is on the extraction of physical quantities without 
simplifying or model assumptions. 
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Figure 9: Experimental data for helicity amplitudes in 
D^K*lv, from [zHH^- 

Units on the vertical axis are 
arbitrary (approximately normalized to unity at f = 0). 



anomalous thresholds are not expected to be large, 
and can be investigated on a case-by-case basis. A 
more complete discussion is beyond the scope of this 



talk and we ignore such complications here. 

In heavy-to-heavy decays it is possible to relate 
form factors by heavy-quark symmetry, making the 
pseudoscalar-vector analysis not significantly differ- 
ent, from the point of view of the number of inde- 
pendent invariant form factors, from the pseudoscalar- 
pseudoscalar case. For heavy-to-light decays, this sim- 
plification is not possible; however, a new symmetry 
emerges due to the large energy of the vector meson, 
and the resulting suppression of the helicity-flip am- 
plitude 0j ■ It is again important to take ad- 
vantage of as much first-principles knowledge as pos- 
sible. By the same variable transformation (jSJ, the 
invariant form factors may all be written in terms of a 
convergent expansion in a small parameter governed 
by the degree of isolation of the semileptonic region. 
Table IIIII shows the maximum size of the parameter 
when to = t+ (1 — ■\/l — t_/t+) , for processes involving 
ground state heavy-light pseudoscalar mesons decay- 
ing into the lowest-lying vector mesons. 

The decay rate can be decomposed in terms of he- 
licity amplitudes 



20 



VtHo _ {l + mL?{EL-rnL')A^-2klA^ 



ruH 



mL(l + 
l + rhL 



(31) 



where hatted variables are in units of mH^ and 
and fci denote the energy and momentum of the light 
vector meson in the rest frame of the heavy pseu- 
doscalar meson. The amplitudes V^Hq, V and Ai 
correspond to cross-channel production of states with 
quantum numbers = 1"*", 1~ and 1"*", respectively. 
On the premise of analyticity, these functions can be 
expanded as in 10}. For the choice of (/>, we use the 
form motivated by unitarity (although the evaluation 
of the unitarity bound itself is not directly relevant). 
Using the results of 9], at Q = 0, 



/z(<,0) 



(t+-to)i/4 V -t 

t- - t 
t+~t f z{t,0) 



to-t 



5/2 



z{t, to) 



-1/2 



t_ - t 



4)Aiit,to) cx 



y/t+ - t f z{t,0)\ f z(t,ta) 



tn-t 



-1/2 



'For the case of _B D*£u, see e.g. 



■'Form factor conventions are as in 143, l7] 
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t- 



t 



-1/4 



(32) 



For the function P, we use for the resonances lying be- 
low threshold: [zllzl = 1" : m = 2.11, 2.63 GeV; 
JP ^1+ ■.m = 2.46, 2.56 GeV. 

Fig. 1^1 shows the resulting invariant form factors ex- 
tracted from the nonparametric analyses in [tqI l80l | , 
after extracting the kinematic function Pcj), and trans- 
forming to the variable z. The small size of z implies 
that these functions should deviate from a straight 
horizontal line by only a few percent over the entire 
range. With this in mind, it is straightforward to 
extract useful information in a systematic way. For 
instance, we can test the ©(Aqcd/^t-c) suppression of 
helicity amplitudes at i = 0: [equivalently, measure 
the form factor ratio y(0)/Ai(0)] 



0.27±0.06 lie 
0.37 ±0.04 E 



(33) 



up to corrections of 0{z). 



8. Summary 

There is interesting information lurking in the 
semileptonic data, with applications to precision phe- 
nomenology, and for exploring a poorly-understood 
limit of QCD. As a practical matter, a simple and 
well-known variable transformation provides a pow- 
erful tool for analyzing and making full use of the 
semileptonic data. The effectiveness of this transfor- 
mation has been obscured by a reliance on unitarity 
bounds, which were used as safe estimates in the ab- 
sence of experimental data. Rigorous arguments for 
a stronger convergence of the expansion |0J) are now 
available, and can be further tested and refined us- 
ing experimental data from many semileptonic decay 
modes. 
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